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ABSTRACT 

Velocity, chemical abundance and spatial distribution data for some 1500 K- and M- 
giants in the Galactic Bulge, which are presented in an accompanying paper (Ibata & 
Gilmore 1995), are analysed. We provide three levels of analysis: an almost model- 
independent determination of basic symmetry properties of the galactic bulge; a more 
detailed analysis based on conservative assumptions, which investigates the basic kine¬ 
matic and chemical abundance distributions in terms of a few parameters; and a more 
detailed investigation of possible functional forms for the density and kinematic dis¬ 
tributions that are consistent with the data using maximum likelihood hts. 

The bulge has a well-determined linear rotation curve over the range TOOpc < 
R < 3500pc with amplitude ~ 25 km s“^ kpc“^. Several velocity dispersion models 
are found to ht the data. We test the oblate isotropic bulge model of Kent (1992) 
by integrating the kinematics predicted by that model into our Galaxy model, and 
comparing with our data. We hnd reasonable agreement except in the most distant 
of our helds from the Galactic centre. We do not hnd a signihcant requirement for 
asymmetry or a bar from our kinematics. 

Our data allow a quantitative comparison between the Milky Way bulge and 
the spiral galaxy bulges studied by Kormendy & Illingworth (1982) over a similar 
Galactocentric distance range; the Galactic bulge appears to be characteristic of the 
extragalactic bulge population. We determine the metallicity distribution of K giants 
in that subset of our helds where the photometric calibration is adequate. Our metal¬ 
licity index measures the Mg& feature, and is calibrated against local held standards. 
Thus our abundance scale assumes a similar relationship between [Mg/Fe] and [Fe/H] 
in the bulge to that which is seen in the solar neighbourhood. Given that assumption, 
the mean metallicities in these helds is [Fe/H] « —0.3 and does not vary signihcantly 
over the regions investigated. This mean is very close to the mean abundance of 
K giants observed recently in Baade’s Window (McWilliam & Rich 1994). That is, 
there is no detectable abundance gradient in the Galactic bulge over the galactocentric 
range 500pc < Rqc < 3500pc. The overall bulge abundance distribution is in good 
agreement with a model that combines standard solar neighbourhood determinations 
of the abundance distribution functions for the halo and thick disk, together with a 
closed-box ‘simple model’ of chemical evolution with yield ~ Q.IZq to describe the 
bulge component. 

We derive the distribution function of specihc angular momentum for the bulge 
from our data, and compare it with determinations for the halo, the thick disk and 
the thin disk from Wyse & Gilmore (1992). We conhrm that the bulge and the halo 
have angular momentum distributions which are indistinguishable, as do the thick 
disk and the thin disk. The bulge-halo distribution is however very different from the 
thick disk-thin disk distribution. This is perhaps the strongest available clue to the 
evolutionary relationships between different Galactic structural components. 

Key words: The Galaxy, Galactic bulge, stellar kinematics, stellar abundances, 
galaxy evolution, galactic structure 


1 THE SPATIAL, KINEMATIC AND 
ABUNDANCE STRUCTURE OF THE BULGE 

The Galactic bulge is a very centrally condensed Galactic 
structural component. Its half-light ‘radius’ along the mi¬ 
nor axis is estimated variously between Si 200 pc for metal 


rich M-giants (Frogel et al.l990) to Si 800 pc for K-band 
surface brightness (Kent et al.l991). The most recent de¬ 
termination, from the DIRBE experiment on GOBE, pro¬ 
vides a scale height of some 300pc (Weiland etal 1994). 
For comparison, the scale height of the Galactic old disk 
is ~ 300pc, while Baade’s Window, the innermost region 
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where the bulge may be well-studied with optical data, is 
some 550pc from the Plane. Power-law hts to inte¬ 
grated surface brightness and to tracer object number den¬ 
sity both require 2.6 < a < 3.7 (Frogel et al.l990, Kent 
et al.l991), with a fa 1.8 inside |6| < 4°. Once the disk 
contribution is deducted, the integrated K-band luminos¬ 
ity (Kent et al.l991), the IRAS point-source surface density 
(Harmon & Gilmore 1988, Weinberg 1992) and the DIRBE 
data display a flattening of c/a fa 0.6. This flattening is 
also seen in the inner RR Lyrae system (Wesselink 1987), 
though the relationship if any between these very old and 
rather metal-poor stars and the dominant population in the 
central bulge, which is rather metal-rich, remains obscure. 

Several lines of evidence show that the bulge is at least 
partly non-axisymmetric: Binney et al.(1991) demonstrated 
that the existence of a bar, extending out to Sf 2.4 kpc from 
the Galactic centre, explains naturally the observed kine¬ 
matics of HI, GO and GS gas in the region \£\ < 10°; fur¬ 
thermore, Blitz & Spergel (1991) showed that the 2.4jUm 
balloon data of Matsumoto et al.(1982) within |6| < 7° is 
brighter at positive £ than at negative £, as would be ex¬ 
pected from a barred stellar distribution. The DIRBE data 
also marginally favour such an asymmetry, which was sus¬ 
pected in the IRAS data (Harmon & Gilmore 1988). From 
an analysis of IRAS point sources selected from those sources 
with |6| < 3°, Weinberg (1992) hnds evidence for a large 
scale stellar bar in the Galactic disk which he suggests ex¬ 
tends out to Sf 5 kpc; however this need not be related to 
the other non-axisymmetric phenomena mentioned above. 

Kinematic studies of candidate bulge tracers (Miras, 
planetary nebulae, K giants, M giants and carbon stars) 
show that the bulge has a central dispersion of Sf 130kms~^ 
which decreases outwards (see e.g. the minor axis data com¬ 
piled by Tyson 1992). In a given region there is a tendency 
for lower abundance tracers to have a higher velocity disper¬ 
sion, as predicted by dissipational collapse models, such as 
the well known ‘ELS’ model (Eggen et al.l962). However, 
available information on the mean rotation of the bulge 
as a function of cylindrical coordinates (R, z) in the Galaxy 
is confusing, because different tracers, or similar tracers with 
different abundance or dispersion, display markedly differ¬ 
ent Vcf,', a compilation of rotation data showing this effect is 
displayed in Figure 1 of Menzies (1990). 

The shape of the bulge velocity ellipsoid has been con¬ 
strained only in Baade’s Window, very near the centre. 
Spaenhauer et al.(1992) obtained proper motions for stars 
in that region; their hndings are consistent with the nuclear 
bulge having nearly isotropic velocity dispersion. 

A recent high resolution spectral survey of bulge K- 
giants in Baade’s Window (McWilliam & Rich 1994) showed 
that the bulge has mean metallicity [Fe/H] = —0.25. That 
work contradicts earlier hndings (e.g. Rich 1988, Frogel & 
Whitford 1987, Geisler & Friel 1992, Tyson 1991) which 
had deduced that the Baade’s Window bulge was extremely 
metal-rich, with a mean abundance [Fe/H] Si -|-0.30. The ex¬ 
istence of an abundance gradient in the inner Galaxy is cur¬ 
rently poorly constrained partly because there is now reason 
to believe that the Washington photometric system, upon 
which many of the relevant studies have been based (Geisler 
& Friel 1992, Tyson 1991, Harding & Morrison 1993), suf¬ 
fers from calibration limitations in the relevant abundance 
range. Nonetheless, there is suggestive evidence for an abun¬ 


dance gradient in the central bulge stellar population in at 
least some tracers: e.g. infrared photometry of M giants is 
consistent with an abundance change of Si 0.4 dex from —3° 
to —12° (Frogel et al.l990). 

1.1 Plan of this Paper 

Our aim in this paper is to present an analysis of our kine¬ 
matic, abundance and photometric data at three comple¬ 
mentary levels. First, in a model - independent way, to 
identify general structural properties of the Galaxy which 
we may be conhdent are reliably determined. Second, in 
a model-dependent way, with a small number of conser¬ 
vative assumptions, to provide a larger amount of detail 
while retaining an appreciation of the importance of the 
various assumptions required. Third, with increased model- 
dependence, to investigate the full amount of information 
potentially contained in the data. This provides a wealth of 
information, though at the expense of risking one’s compre¬ 
hension in a sea of parameters. 

In §2 we briefly summarise the data which is presented 
in detail in an accompanying paper (Ibata & Gilmore 1995). 
§3 outlines the stellar luminosity classihcation and describes 
the implications of the reliability of this classihcation for 
the astrophysical conclusions of this paper. §4 presents the 
most robust and least model-dependent deductions from the 
data, concluding that the bulge is kinematically symmetric, 
in spite of the confusing presence of the Sagittarius dwarf 
galaxy in our data. §5 discusses our photometric star count 
data, together with other surface brightness measures, to 
deduce the most appropriate bulge and old disk spatial den¬ 
sity distributions. Surprisingly, our data suggest a low inner 
disk scale height. §6 provides a kinematic analysis utilis¬ 
ing Kent’s isotropic oblate rotator dynamical model and the 
density distribution of §5. The agreement of our data with 
the predictions of this model, at least in the inner few kpc 
of the Galaxy, is rather good. In §7 we extend the kinemat- 
ical analysis to a more wide-ranging investigation of possi¬ 
ble kinematic models, rotation curves, dispersion prohles, 
and velocity ellipsoid orientations and shapes. We employ 
a maximum likelihood method to compare models with the 
data, and analyse the data both held-by-held and en masse. 
A robust conclusion is that the bulge has a linear rotation 
curve with well-determined amplitude over the whole range 
of our data. §8 takes up this robust determination, and 
derives the corresponding distribution function of specihc 
angular momentum for the bulge. We compare this to simi¬ 
lar distribution functions for other Galactic components, the 
halo, the thick disk and the thin disk, to show how the bulge 
is closely related to the halo, but not to the disks. In §9 we 
combine the abundance and kinematic data, and use these 
to isolate the abundance distribution of the bulge. We then 
compare this to the ‘closed-box simple model’ of chemical 
evolution, generating a rather good ht with plausible param¬ 
eters. In §10 we compare our bulge with data for external 
galaxies, which we are able to do over a similar range of 
galacto-centric distances for the hrst time. §11 summarises 
and concludes. 

2 SUMMARY OF THE DATA 

In an accompanying paper, Ibata & Gilmore (1995, hence¬ 
forth Paper 1) present the results of an analysis of photo- 
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metric and spectroscopic data for a large sample of bulge 
giants. Tire stellar data are presented in Table B.l of Paper 
f in detail. Each star has a position, a colour, a velocity, 
a metallicity estimate, a spectral classihcation and an es¬ 
timate of its luminosity class (dwarf or giant). A total of 
Si 1500 stars were observed spectroscopically. In summary: 

i) Five carbon stars were discovered in the spectroscopic 
survey. Four of these stars lie in the i = 5° helds, are 
members of the Sagittarius dwarf galaxy, and are discussed 
further in Ibata, Gilmore & Irwin (1995; hereafter IGl). 

ii) Forty-two M giants were found in the complete sample, 
of which 36 also lie in the 1 = 5° helds. Again, they are 
members of the Sagittarius dwarf, and are discussed further 
by IGl. 

hi) Some 250 K dwarfs (si 20% of the sample) were classi- 
hed and identihed both by visual inspection of the spectra, 
and by use of a new classihcation technique based on Prin¬ 
cipal Gomponent Analysis. This technique will be discussed 
further elesewhere (Ibata & Irwin, in preparation). 

iv) The stars of direct interest as members of the Galactic 
bulge, the K giants, numbered Si 1200, making up Si 80% of 
the observed sample. 

3 DWARF CONTAMINATION 

Since we wish to determine parametric descriptions of the 
kinematics of the Galactic bulge, it is important that any 
stars in our sample which lie in the foreground disk or the 
background Sagittarius dwarf galaxy are reliably identihed. 
For present purposes, we describe such stars as ‘contami¬ 
nants’. Dwarf-giant luminosity classihcation has been car¬ 
ried out for all our stars. Giant carbon stars and M-giants 
can be classihed very reliably (cf Paper 1), in large part since 
they are apparently bright, and so have high signal-noise 
ratio spectra. The fainter K dwarfs may however be less 
reliably classihed. Thus we check by testing for their ehect 
on the abundance and kinematic distributions. Before this, 
however, we emphasise, as discussed in paper 1, that the 
number of stars in our sample classihed as disk K dwarfs 
is in excellent agreement with the number predicted from 
our Galaxy model. There is therefore no reason to expect 
that any signihcant error in the luminosity classihcation has 
occured. 

3.1 Contamination in the Velocity Distribntions 

The reliability of K dwarf - giant discrimination signihcantly 
ahects the conclusions that will be drawn from this data set. 
We therefore hrst investigate whether the objects identihed 
as K dwarfs have an expected velocity signature. Any K 
dwarfs in the sample will be foreground disk stars, some 
Ikpc from the Sun (cf Paper 1). They will therefore have 
a clear kinematic distribution characteristic of the old disk. 
The observed velocity distributions of all stars classihed as 
dwarfs are plotted in Figure 1, where we have superimposed 
the expected distribution according to our galaxy model. All 
hts are good to better than the la level, so we are conhdent 
that the classihcation scheme does indeed pick out dwarfs, 
and does not select a signihcant number of giants. The 
relationship of disk and bulge kinematics down these lines 
of sight, which determines the ehect of any possible residual 
dwarf contamination in the giant sample, is discussed later. 


3.2 Contamination in the Metallicity 
Distribntions 

The stellar sample was divided into into giants and dwarfs 
using a Principal Gomponent Analysis technique (Ibata & 
Irwin in preparation). This classihcation procedure is a 
function of the colour of the star and a coefhcient similar 
to its linestrength The metallicity estimates we obtained (cf 
Paper 1) are the result of an interpolation in an empirical 
correlation between colour and linestrength, calibrated by 
data for many standard stars. Thus we expect our classi¬ 
hcation and metallicity estimates to be strongly correlated. 
We note that the ‘metallicity’ derived for a dwarf is not cal¬ 
ibrated in this project, so that such ‘metallicity’ estimates 
do not correspond to meaningful abundances. In particular, 
no zero point calibration of the dwarf abundance estimates 
has been attempted. 

We examine the shape of the (pseudo-)‘metallicity’ dis¬ 
tribution of stars classihed as dwarfs, so as to see the way 
that the K giant metallicity distributions would be con¬ 
taminated by any residual mis-classihcation. The distribu¬ 
tions are shown in Figure 2; objects classihed as dwarfs are 
strongly concentrated in the region [Fe/H] > 0. Thus, any 
errors in luminosity classihcation will preferentially provide 
contamination to the K giant metallicity distribution also in 
that range. Recall also (paper 1) that the metallicity distri¬ 
butions are reliable only in those regions where a reliable cal¬ 
ibration to the APM photometry is available — that is, not 
in the helds at (£ = —15°, 6 = —12°) or (£ = 5°, 6 = —20°), 
( cf Paper 1). 

4 MODEL-INDEPENDENT DIFFERENCES IN 
THE VELOCITY DISTRIBUTIONS 

Using the non-parametric Kolmogorov-Smirnov (KS) test 
we compare the observed radial velocity distributions of dif¬ 
ferent parts of the Galactic bulge. A cumulative probability 
distribution is constructed from the radial velocity data set 
in each region such that the probability of the lower veloc¬ 
ity bound is zero and the probability of the upper bound is 
unity (the value of these bounds will be stated later). The 
KS statistic D is found: this is simply the maximum value 
of the absolute diherence between the two cumulative prob¬ 
ability distributions. The KS test allows one to calculate 
the probability P(D > observed) that D should be greater 
than the observed value if the distributions are drawn ran¬ 
domly from the same parent population. If this probability 
is either very small or very near unity then the distributions 
are signihcantly different. 

The hrst question to be addressed is: is there any ev¬ 
idence that the observed bulge-region stars are kinemati¬ 
cally non-axisymmetric? This question can be answered di¬ 
rectly (in a model independent way) by comparing pairs of 
lines of sight which are situated at points rehected symmet¬ 
rically about the bulge minor axis. The regions that we 
can test in this way are those at (£ = —5°, 6 = —12°) and 
(£ = -|-5°, h = —12°). Then the above question can be stated 
as: if the map v —v is applied to the velocity distribu¬ 
tions at negative longitudes, is P(D > observed) such that 
the distributions are signihcantly diherent? 

The answer to this depends on the choice of velocity 
bounds for the KS test. Two distributions are formed by 
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Figure 1. Comparison between the expected K dwarf velocity distributions (smooth curves) and those for the stars classified as K 
dwarfs. The agreement between the data and the model prediction lends confidence to the classification. The lines of sight are: (a) 

e = -25°, b = -12°, (b) e = -15°, b = -12°, (c) e = - 5 °, b = - 12 °, (d) e = 5°, b = - 12 °, (e) e = 5°, b = - 15 ° and (f) e = 5°, 
b = -20°. 


selecting stars from the {£ = —5°, 6 = —12°) and {£ = 
+5°, 6 = —12°) fields starting from a lower velocity bound 
of —250kms~^ and ending at an upper velocity bound v. 
We then calculate the probability P{D > Observed), and 
equate this probability with the probability that an event 
should be found more than x standard deviations from the 


mean of a gaussian distribution of standard deviation a: 

+ CO 

2 _ / exp( — t^/2(7^)(it = P(/? > Observed). (1) 

v27r J 

X 

Thus solving for x, we find the probability P{D > Observed) 
in units of cr; this is shown in Figure 3. Clearly, up to 
V Si 170kms“^ the radial velocity distributions in the two 
regions are distinguishable only below the 2a level, which 



The Outer Galactic Bulge: II - Analysis 5 


Figure 2. The pseudo-‘metallicity’ distribution of those stars classified as dwarfs is shown for (a) i = —25°, b = —12°, (b) f = —15°, 
b = —12°, (c) I = —5°, b = —12°, (d) I = 5°, b = —12°, (e) I = 5°, b = —15° and (f) I = 5°, b = —20°. These abundances are neither 
calibrated nor reliable. Their significance is to show that any residual dwarf contamination in the K giant star sample will have greatest 
effect at the high metallicity end of the K-giant distribution, for [Fe/H] ^ 0.0. 


is not significant, but beyond v Si 170kms~^ they are sig¬ 
nificantly (si 4.5a) different. The (statistically significant!) 
kinematic feature at Si 170kms~^ in the 5° fields is due to 
the (hitherto unknown) Sagittarius dwarf galaxy evident in 
the £ = -|-5° fields (Ibata, Gilmore & Irwin 1994; IGI). It 
is the presence of the Sagittarius dwarf in our survey fields 
which leads to the need for the relatively complex consid¬ 
eration of velocity bounds in Figure 3. At velocities be¬ 
low ~ 170km/s, where the data exclusively represent the 


Galactic bulge, the two symmetric fields do not differ sig¬ 
nificantly in their kinematics (Figure 3). Thus we conclude 
that, withing the precision of the present data, the bulge is 
kinematically symmetric. 

The next question to be asked is: is there evidence for 
a change in kinematics with (vertical) distance from the 
Galactic plane in the regions investigated? This may be 
tested directly with the radial velocity data obtained from 
the three fields parallel to the Galactic minor axis, {£ = 
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Figure 3. Testing for any possible difference between the 
(e = —5°,b = —12°) and (£ = +5°, 6 = —12°) fields. Two cu¬ 
mulative distributions are formed between —250kms“^ and an 
upper bound v from the two data sets. The KS test statistic D, 
(the maximum value of the absolute difference between the two 
cumulative distributions) is found. The plot shows the probabil¬ 
ity P[D > Observed) (in units of c — see the text) as a function 
of the upper bound v. Note that below v % 170kms“^, the two 
distributions cannot be considered different. The velocities shown 
are ‘Galactocentric radial velocities’ — this term is explained in 
Paper I. 

5°,6 = -12°), {£ = 5°,6 = -15°) and {£ = 5°,6 = -20°). 
It is found that, constrained by any upper or lower velocity 
bounds, the distributions are different at less than 1.5(7. Fit¬ 
ting a gaussian to the data in the heliocentric radial velocity 
range —lOOkms”^ <v< -flOOkms”^ (so as to avoid bias¬ 
ing the ht with stars from the Sagittarius dwarf galaxy), we 
hnd that the mean heliocentric velocities of the maximum- 
likelihood gaussian hts in the helds at 6 = —12°, h = —15° 
and b = —20° are —6.8 ± 4.1kms“^, —6.4 ± 5.6kms“^ 
and —10.7 ± 7.0kms“^ respectively. We therefore conclude 
that there is no evidence from the data for vertical varia¬ 
tions in the kinematics of bulge-region stars in the range 
— 1.6 ^ z > —2.7kpc at £ = 5°. Note that this does not 
necessarily imply that there is no vertical variation in bulge 
kinematics, since we have not considered the expected kine¬ 
matics of the disk, thick disk and halo populations yet. 

It has been noted that bulges of external galaxies that 
have so called ‘boxy’ isophotes generally display cylindrical 
rotation (Kormendy & Illingworth 1982, Shaw et al.l993); 
this phenomenon is observed in many low luminosity bulges 
(Jarvis 1986). Kent et al.(1991) hnds that the 2.2jum Galac¬ 
tic bulge emission has slightly boxy contours, which has sug¬ 
gested the possibility of cylindrical rotation. The recent 
DIRBE data however (Weiland etal 1994) show this shape 
to be an artefact of patchy reddening. Since both the mean 
velocity and the velocity distributions in our helds at 5° ap¬ 
pear to be indistinguishable, over a factor of two in distance 
from the plane, our data are consistent with the bulge being 
a cylindrical rotator. 

It should be noted that the results presented in this sec¬ 
tion should be considered only as suggestive. This is because 
we have not taken the selection function into account in our 
calculations: stars in the different helds were not selected 
in an identical way. Thus even if the kinematics of stars in 
two helds were intrinsically identical, the diherences in the 


selection functions would give rise to diherent proportions of 
disk, thick disk and halo stars in the two helds, and one ob¬ 
tains a spurious comparison. In section 7 we shall model the 
kinematic distributions and amend this dehciency by taking 
the selection function into account. Since this sophistica¬ 
tion in modelling will come at the cost of introducing many 
free parameters into the discussion, we believe this simplistic 
discussion valuable to complement the later complexity. 

5 BULGE MODEL STAR COUNTS 

The contributions to the total kinematics from the thin disk, 
thick disk and halo are kept hxed (with the exception of the 
old disk scale height) as given by the parameters and rela¬ 
tions detailed in our star count model of the Galaxy. This 
model will be described more fully in its present form else¬ 
where (Gilmore, Wyse & Hewett in preparation), but is ba¬ 
sically that described by Gilmore (1984), by Gilmore, Wyse 

6 Kuijken (1989), and by Gilmore, King & van der Kruit 
(1990). The important new features here are the parame¬ 
ters specihc to the Galactic bulge, which we outline in this 
section. 

The luminosity function (LF) adopted for the bulge 
component is taken from the disk LF by Wielen (1974) with 
the LF for M-giants in Baade’s Window by Frogel & Whit- 
ford (1987) appended on at the bright end. An arbitrary 
normalisation is used, because the zero-point of the bulge 
density law is not known. 

The colour-magnitude (GM) relation used to represent 
the bulge is that of M67 compiled by Ghiu (1980). This disk 
cluster is of slightly higher metallicity ([Fe/H] = —0.04) than 
the bulk of our K-giant sample (see section 9). 

5.1 Star Count Data and the Galaxy Model 

We investigate the ability of the Galaxy model to pre¬ 
dict stellar number counts given the selection function im¬ 
posed on our data. Stars were chosen for observation 
(c/ Paper I) from the colour-magnitude parameter space 
14.5 < R < 16.5, Bj — R > 1.4; the expected Bj — R distri¬ 
bution of these stars in the {£ = —5°, 6 = —12°) held is 
shown in Figure 4, where we have superimposed the ob¬ 
served distribution. In this simulation, the model works 
acceptably well in the red (beyond Bj — R = 1.4, where our 
data lie), but gives a poor ht in the blue. The steep gradient 
on the red side of Figure 4 clearly shows that it is necessary 
to calibrate the zero-point of the Bj, R photometry to better 
than Si 0.1™ in order to obtain star counts accurate to bet¬ 
ter than Si 30%. Recall that in the (£ = —15°, 6 = —12°), 
(£ = —5°, b = —12°) and (£ = 5° ,b = —20°) helds no reliable 
GGD calibration is available. Reddening errors of Si 0.1™ 
in the Burstein & Heiles maps will introduce similar Si 30% 
errors in the starcounts.Since we have neither precise photo¬ 
metric data nor accurate reddening determinations in most 
of our helds, the data presented here are not well suited for 
the study of stellar number density. 

5.2 Bulge Density Laws 

Given the assumed luminosity function and colour magni¬ 
tude relation for the bulge, a density relation n must be 
chosen to ht the starcounts. However, we show in IGI 
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Figure 4. The Bj — R distribution of stars along the line of sight 
(£ = —5°,6 = —12°) in the magnitude interval 14.5 <R< 16.5 is 
shown: the ‘star’ symbols represent the expected distribution, 
while ‘dots’ represent the observed distribution. 


that the stars with heliocentric radial velocities in the range 
V = 140 ± 10 km s“^ in the 7 = 5° helds belong to the Sagit¬ 
tarius dwarf galaxy, a population that is not accounted for 
in the galaxy model; in those helds the ht of the model 
to the starcounts will therefore be performed on stars with 
heliocentric radial velocities v < 120kms~^. With this re¬ 
striction, we calculate the total number of stars expected 
M, and observed N, in each held i. The zero point of the 
adopted bulge density relation is obtained by minimising 
-N,f. 

Kent et al (1991) modeled 2.4jUm SPAC'ELAB data 
(Kent et al 1992); in the regions relevant to this investiga¬ 
tion, the ‘best ht’ density relation they hnd is p oc /To(s/667) 
(this gives boxy projected density contours), where Ko is a 
modihed Bessel function and s = (R'* + (^/O.bl)’')^^’'. The 
‘best ht’ also contains a double exponential disk with a scale 
length tiR = 3001 pc and a scale height which is a constant 
hz = 165 pc from R = 0 to R = 5300 pc but beyond which it 
increases linearly so as to be hz = 247 pc at Rq. The prin¬ 
cipal contributors to the 2.4jUm emission are late K- and 
M-giants (Jones et al 1981), so we convert standard K-band 
hux into K- and M-giant number counts using the K-band 
luminosity function by Jones et al (1981). The combined 
model for both bulge and disk works acceptably well (c/ 
section 5.1); we hnd that the expected number counts do 
not diher from the observed number counts by more than 
35% in any of the helds. This conhrms that the data of Kent 
et al (1992) is not dominated by a few very bright M-giants. 

Sellwood & Sanders (1988) ht a ‘maximum disk’ model 
to 2.4jUm balloon data by Matsumoto et al (1982) and to 
2.2jum data by Becklin & Neugebauer (1968). The bulge 
density law they hnd is p oc r~^ ®[l + (r/lkpc]~^], which 
clearly asymptotes to We therefore also try an 

like model; then the best ht to K-giant number density has 
(s/so)”^'^, where = R^ + {z/0.6)^ and so = 900pc. 

However, we hnd that, as long as the old disk contribu¬ 
tion to the total K-giant number counts is less than Si 30% 
(zh < 220pc), diherent bulge density laws — that are able 
to ht the number counts — do not signihcantly ahect the 
bulge kinematic parameters derived below. 


5.3 The Old Disk Scale Height 

The disk components of spiral galaxies show constant verti¬ 
cal scale height hz independent of R (see e.g. the chapters 
by van der Kruit in Gilmore, King & van der Kruit 1990). 
Assuming that the presence of the bulge has no ehect on the 
disk scale height in the central regions of the Milky Way, one 
would expect Zh = constant through the Milky Way. Then 
Zh can be set by the solar neighbourhood value: for instance 
Kuijken & Gilmore 1989 ht hz = 249 pc to local K-dwarfs. 
Kent et al (1991) on the other hand hnd a best ht to 2.4jUm 
photometry with a bulge and a disk model in which the disk 
scale height is constant (zh = 165 pc) out to a certain radius 
(R = 5300pc), but beyond which it increases linearly out¬ 
wards so as to be hz = 250 pc at the solar neighbourhood. 
The obvious interpretation is that the 2.4jUm data are dom¬ 
inated by the young disk, which has a scale height < 150pc, 
and not by the old disk. We now check to see if a similar 
ehect appears in our data. 

In the regions investigated, the ratio of expected num¬ 
ber counts between a model with hz = 165 pc and one with 
hz = 300 pc is Si 0.13; this makes the diherence between 
a small disk contribution and one that has, in some helds, 
more than twice the observed number of stars. Since we have 
three helds at 7 = 5°, we can easily compare the vertical be¬ 
haviour of the disk density at R = 700 pc. We deduce that, 
if Zh ^ 300 pc, then either the adopted Wielen (1974) lumi¬ 
nosity function or the colour-magnitude relation of M67 are 
inappropriate approximations for the inner parts of the disk, 
or there are signihcant zero-point diherences in the adopted 
photometry (c/Paper 1) between the helds. The problems 
with such a large scale height are demonstrated in Figure 
5. For 165 < Zh ^ 250 pc, it is possible to ht the 7 = 5° 
helds to within the Si 30% starcount error (c/section 5.1); 
but A z P 220 pc, we are unable to ht the kinematics of the 
held at (7 = —25°, 6 = —12°) at better that about the do- 
level, unless the bulge velocity dispersion in that held has 
the unrealistically low value of Si 30kms~^, when the ht 
is acceptable at the 2o- level. However, these problems are 
resolved if we adopt a scale height Zh ^ 200 pc, and we can 
make acceptable hts to both the starcounts and kinematics 
in all helds. 

Thus in this respect, it appears that the disk of the 
Milky Way is not representative of the disk populations of 
spiral galaxies. The extent to which this is a structural 
feature of the Galaxy, indicates that our survey is still over¬ 
sensitive to young disk rather than old disk stars, or indi¬ 
cates limitations in the calibration of our photometric data 
remains unclear. 


6 COMPARISON TO KENT’S ‘ISOTROPIC 
OBLATE ROTATOR’ BULGE MODEL 

The only kinematic model to date that has been constructed 
from a ht to the large-scale mass structure of the Galactic 
bulge and inner Galactic disk (as obtained from 2.4jUm emis¬ 
sion with an assumption of the mass to light ratio) is that 
of Kent (1992). He assumes that the bulge is an axisym- 
metric, oblate spheroid, of constant M/L and with isotropic 
velocity dispersion, so that the distribution function / is a 
function only of total energy E and the 2 —component of 
angular momentum Lz- The Jeans’ equations then take the 
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Figure 5. Comparison between the observed velocity distributions and the Galaxy model, where we have adopted a scale height of 
Zfi = 300pc for the disk component. The lines of sight are: (a) i = —25°, h = —12°, (b) i = —15°, h = —12°, (c) i = —5°, h = —12°, 
(d) I = 5°, b = —12°, (e) i = 5°, b = —15° and (f) i = 5°, b = —20°. In each panel, the dashed-single-dotted curve represents the halo, 
the dotted curve represents the thick disk, the dashed-triple-dotted curve represents the bulge and the dashed curve represents the disk. 
We have normalised the counts so as to fit the distribution in (f). The fits in (a) and (b) can be rejected at about the 5 (t level, while (c) 
and (d) have too many stars by CS 70% and can be rejected at the Scr level. 


form: 


(TRiRz) = al 



9 ^- , 


o S'!/ 


R djnaR) 
n dR 


where (R, (j), z) are cylindrical coordinates, n is stellar num¬ 
ber density, dl is the Galactic potential, is the mean ro- 
(2a) tational velocity of the bulge and a is the bulge velocity 

dispersion. 

The forces ST /dR and ST /dz are computed, using a 
technique described in Kent (1989), from the mass distribu- 
(2fc) tion of the disk and bulge, as found by Kent et al.(1991). 
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with the assumption that M/L = 1. A central 10® Mq black 
hole is also included in the model, though is not relevant for 
the present study. An algorithm to calculate Sd/ /dR and 
Sd /Idz was kindly made available to us by Steve Kent (pri¬ 
vate communication). The mean velocity and dispersion 
resulting from equation 2 are in good agreement with obser¬ 
vations for \l\ A 100pc and |6| < Ikpc, ie, interior to almost 
all our data, but has hitherto not been tested in more distant 
regions from the Galactic centre. 

The velocity distribution of the bulge population at 
each point along the line of sight is contructed from the esti¬ 
mates of a{R,z) and obtained from equations 2. We use 
the Galaxy model to integrate disk, thick disk and halo con¬ 
tributions, and add Kent’s model for the bulge component; 
this gives the distributions shown in Figure 6. According to 
KS tests, all the hts except that in Figure 6e are unaccept¬ 
able at about the 2-3(7 level. 

Since the observed distributions are more peaked than 
the above model predicts, we ask the question: are these 
bad hts due to contamination from the local dwarf popula¬ 
tion? Objects classihed as dwarfs contaminate the sample 
mostly at the high metallicity end, as we showed in Figure 2. 
To investigate this question, we apply the above kinematic 
model to only those stars with metallicities [Fe/H] < —0.5. 
The resulting hts are shown in Figure 7; now all are accept¬ 
able at better than the 1.5(7 level, except for the held at 
(£ = —25°, 6 = —12°), which can be discarded at the fa 4.5(7 
level. Thus Kent’s model does not work well at Galactocen- 
tric distances R fa 3.5 kpc. It also appears that either the 
model does not predict sufhcient numbers of dwarf stars, 
the apparent excellent agreement between model prediction 
and classihcation of dwarfs noted in §3 is chance and that 
there is signhcant residual dwarf contamination in the giant 
sample, or that there is a signihcant dependence of velocity 
dispersion on metallicity in the bulge population — these 
possibilities will be investigated more thoroughly in section 
9 below. 


7 FUNCTIONAL FITS TO BULGE 
KINEMATICS 

7.1 Maximum Likelihood 

When examining the models that can be allowed or refuted 
by the data we hx most of the relations in the galaxy model. 
We wish to hnd the most likely parameters of the remaining 
relations. The likelihood L{H\D) of a hypothesis (or model) 
H given data D is dehned to be proportional to the prob¬ 
ability P{D\H) of the data given the hypothesis (see e.g. 
Edwards 1992). Maximising the likelihood L 

n 

L<^\{P(D,\H), (3) 

^=l 

is equivalent to minimising the quantity 
s = -2 ^lnP(U,|W), (4) 

i=l,n 

where D = {Di,D 2 , ■ ■ ■ ,Dn}- (It is customary to include 
the factor 2 in the dehnition of s so as to emphasize the 
analogy with thermodynamic entropy). Given the observed 


{i = 1, . . . , m} helds containing rit stars, the statistic s is 
found from 


s 


-2 


m lit 


(=1 j=i 


Gjvl) 

T ’ 


(5) 


where G is the expected line of sight velocity distribution, 
T, = f G(v)dv and Vj are the velocity values. 

The statistic s is minimised by varying the chosen pa¬ 
rameters in the model with an ‘amoeba’ minimising algo¬ 
rithm from Press et al (1986). 


7.2 Bootstrapping 

With conventional Monte-Garlo simulation, one estimates 
the conhdence limits on the parameters ao htted to a data 
set by Ending new parameters a, which are fitted to a large 
number i of synthetic data sets. The expectation is that the 
probability distribution of a, — ao mirrors a, — atrue, where 
atrue are the parameters one would find from a hypothet¬ 
ical infinitely precise experiment. The synthetic data sets 
are created according to one’s understanding of the process 
being investigated and the related measurement errors. 

However, since it is diflfcult to quantify the errors in the 
input parameters to our starcount model, we use instead the 
so-called bootstrap method, which is a Monte-Garlo simu¬ 
lation where the synthetic data sets are created by drawing 
N points, with replacement, from the original data set of 
N points. Taking i w 50 synthetic data sets in each of the 
observed lines of sight, we sample the distribution a, — atrue 
from the i maximum likelihood fits a,. For a discussion of 
this technique, see Press et al (1986). 


7.3 The Fits 

We now investigate some possible functional forms of cr{R, z) 
and in the bulge, finding parameters that are consistent 
with our data. The most-likely bulge rotation relations given 
the data are dependent on the bulge velocity dispersion re¬ 
lations and upon the scale height zu of the disk component 
in the inner regions of the Galaxy. We argued above (in sec¬ 
tion 5.3) our preference for the values zu ^ 200 pc. Given 
that this is a major source of uncertainty in the fits below, 
we will present three possibilities with each ht: (N) - no 
disk, (K) - for zu = 165 pc (as given by Kent et al.l991) - 
for Zh = 200 pc (the largest scale height that appears to be 
consistent with our data). The KS test probabilities for all 
the fits discussed below are given in table 1; except when 
otherwise stated, they are all acceptable at better than the 
Si 2(7 level. 

Each of the observed velocity distributions has infor¬ 
mation on the kinematics of the Galaxy over a large range 
of Galactocentric radius. Thus the kinematic model of the 
Galaxy should be fit simultaneously to the velocity distribu¬ 
tions in all of our fields. However, the maximum likelihood 
fits will naturally be biased towards finding parameters that 
fit the better sampled fields best. We therefore first in¬ 
vestigate, in the next two paragraphs below, what can be 
said about linear and cylindrical rotation in the bulge, if 
the kinematic distributions from each held are fit indepen¬ 
dently from one another. The results of these hts can later 
be compared to those in which the velocity distributions are 
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Figure 6. Comparison between the observed velocity distributions and the expected distribution according to a model which has 
standard parameters for the disk, thick disk and halo populations, but whose bulge population is an oblate isotropic rotator as fit by 
Kent (1992). Kent’s density model for the bulge and thin disk works reasonably well, giving number counts that agree with observations 
to % 35%. Note that the total counts in the modeled distributions below have been normalised to the counts in the data, and not 
predicted. The lines of sight are: (a) t = —25°, b = —12°, (b) t = —15°, b = —12°, (c) t = —5°, b = —12°, (d) t = 5°, b = —12°, (e) 
I = 5°, b = —15° and (f) 1 = 5°, b = —20°. The fit in (e) is acceptable at the la level, though all the rest are poor, at about the 3a 
level. 


fit simultaneously. 

1. Oblate bulge 

We investigate models in which the bulge has an 
oblate stellar density distribution of the form n{R, z) = 


r'o(s/so)~^'^, where + {z/0.6)^ and so = 900pc 

(c/ section 5.3). The KS test comparisons between kine¬ 
matic distributions from data on opposite sides of the bulge 
strongly support this model (c/section 4). 

l.a. Is the bulge a linear rotator? 
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Figure 7. The velocity distributions of K giant stars with metallicities [Fe/H] < —0.5 are compared to a model similar to that shown 
in Figure 6, but without a local dwarf population. The lines of sight are: (a) I = —25°, h = —12°, (b) I = —15°, h = —12°, (c) I = —5°, 
h = —12°, (d) I = 5°, h = —12°, (e) I = 5°, h = —15° and (f) I = 5°, h = —20°. All fits, except that in (a), are acceptable at better 
than the 1.5(T level. The fit in (a) can be discarded at the % 4.5(T level. This model therefore appears to work well in the central parts 
of the Galaxy. 


The velocity distribution in each of our helds contains 
information about the mean azimuthal velocity and the ve¬ 
locity dispersion tensor over a large range of heliocentric 
distance. However, the largest contribution to the distribu¬ 
tion will come from regions near the tangent point to the 
line of sight in the particular held under study, where the 
stellar number density is highest. In a small distance range 
near the tangent point, we can expect the mean azimuthal 


velocity to be a linear function of Galactocentric radius 
R, so to hrst order we can attempt to ht with a func¬ 
tion v^(R) = HbR. We proceed to ht the coefhcient Qb 
independently for each held. Then, if the htted parame¬ 
ters Qb are the same over all helds (to within the measuring 
errors), the bulge can indeed be said to be a linear rota¬ 
tor. For these hts we select only those stars classihed as 
giants with metallicities [Fe/H] < —0.5; this should allevi- 
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ate contamination problems with local K dwarfs (see Figure 
2) and distant disk giants (whose mean metallicity we can 
expect to be approximately solar). We also select only the 
Helds at (I = -25°, 6 = -12°), (I = -15°, 6 = -12°) and 
{I = -5°, 6 = -12°). The fields at [I = 5°, h = -15°) 
and {£ = 5°, 6 = —20°) are not analysed because we wish 
to suppress vertical terms in the fits, while the field at 
{£ = 5°, 6 = —12°) has a velocity distribution that is sub¬ 
stantially contaminated with stars belonging to the Sagit¬ 
tarius dwarf galaxy (this is discussed in IGl). We will as¬ 
sume that the R component of the velocity dispersion is 
almost constant in the neighbourhood of the tangent point, 
i.e. <trr(R) Si constant. 

l.a.i. Assuming that the bulge has isotropic velocity dis¬ 
persion, i.e. (Trr = = (Tzz, we find 

(1=—25°,6=—12°): Qb = 26.6± 3.1 km s~^kpc~^ and 
(7 = 48 ± 3km s“\ 

(1=—15°,6=—12°): Qb = 25.3± 4.5 km s~^kpc~^ and 
(7 = 45 ± 6km s“\ 

{£ = -5°,6 = -12°): Qb = 24.8 ± 2.7kms“^kpc“^ and 
(7 = 71 ± 4km s“\ 

l.a.ii. Constraining <trr : : Uzz in the ratio of solar 

neighbourhood disk dispersions, 39 : 23 : 20, gives: 

{£ = —25°, h = —12°): Qb = 26.6 ± 3.7km s~^kpc~^ and 
(7 = 74 ± 5 km s~^, 

(£ = -15°, 6 = -12°): Qb = 25.3 ± 4.3kms“^kpc“^ and 
(7 = 67 ± lOkms”^, 

(1=—5°,6=—12°): Qb = 25.4±2.9kms~^kpc~^ and 
(7 = 86 ± 7km s~^, 

l.a.iii. Constraining ctrr : : Uzz in the ratio of solar 

neighbourhood halo dispersions, 131 : 102 : 89, gives: 

{£ = -25°, 6 = -12°): Qb = 26.6 ± 2.9kms“^kpc“^ and 
(7 = 59 ± 4km s“\ 

{£ = —15°, 6 = —12°): Qb = 25.1 ± 5.2kms“^kpc“^ and 
(7 = 61 ± 6km s“\ 

(1=—5°,6=—12°): Qb = 25.1±2.0kms~^kpc~^ and 
(7 = 77 ± 5 km s“^, 

Clearly, Qb is not very sensitive to the assumed velocity 
dispersion model (the detailed relation between Qb and a 
will be explored below). Thus to good approximation, the 
bulge, as sampled with our data set, is a linear rotator, 
whose rotation rate is Qb ~ 25 km s~^kpc~^. It is interesting 
to note that, extrapolated to the solar neighbourhood, this 
population would have a mean azimuthal velocity identical 
to that of the solar neighbourhood thin disk, though this is 
perhaps coincidental. 

l.b. Is the bulge a cylindrical rotator? 

Some calculations presented in section 4 suggested that the 
mean rotation of the bulge does not vary significantly with 
height above the Galactic plane. We now re-investigate this 
point, using data from the three fields at (1 = 5°, 6 = —12°), 
(1 = 5°,6 = —15°) and [£ = 5°,6 = —20°). We would ide¬ 
ally like to choose, as above, only those stars with (say) 
[Fe/H] < —0.5, that way the contribution of the local and 
distant disk populations to any vertical variations in the 
mean velocity would be reduced. Unfortunately, we have 
no reliable photometry and hence no reliable classifications 


or metallicities for stars in the field at {£ = 5°, 6 = —20°); 
however, comparing the {£ = 5°, 6 = —12°) field to that at 
{£ = 5°, 6 = —20°) gives a large vertical baseline, which is 
crucial in answering the cylindrical rotation question. Fur¬ 
thermore, if we accept the preliminary photometric calibra¬ 
tions for the (£= 5°, b = —20°) field (c/Paper 1), and choose 
only those stars with [Fe/H] < —0.5, we would obtain the 
velocity distributions shown in Figure 7 which have too few 
stars for reliable modeling. 

We will fit the velocity distribution (of all the stars) of 
one field at a time, finding fit and a for each field. (In the 
next section we will attempt to fit fit and a simultaneously 
in all our fields). Again we assume that crrIR) is almost 
constant in the neighbourhood of the tangent point. 

l.b.i. Assuming that the velocity ellipsoid is isotropic, i.e. 
IJRR = <^ 4 , 4 , = Uzz, we find 

(N) (£ = 5°, 6 = —12°): Qb = 15 ± 8kms“^kpc“^ and 
(7 = 54 ± 5 km s~^, 

(£= 5°, b = —15°): Qb = 22 ± 12 km s~^kpc~^ and 
(7 = 70 ± 7km s“^, 

(f=5°,6=—20°): Qb = 18±10kms“^kpc“^ and 
(7 = 42 ± 5 km s~^. 

(K) (f = 5°,6 = —12°): Qb = 14±6kms~^kpc~^ and 

(7 = 54 ± 5 km s~^, 

(£= 5°, b = —15°): Qb = 17 ± 11 km s~^kpc~^ and 
(7 = 67 ± 8 km s~^, 

(£= 5°, b = —20°): Qb = —7 ± 15 km s~^kpc~^ and 
(7 = 44 ± 6 km s“^. 

(L) (f = 5°,6 = —12°): Qb = 21±5kms~^kpc~^ and 

(7 = 49 ± 4 km s“^, 

(f=5°,6=—15°): Qb = 16±8kms~^kpc~^ and 
(7 = 66 ± 10 km s“\ 

(f=5°,6=—20°): Qb = —8±10kms“^kpc“^ and 
(7 = 44 ± 9 km s“^. 

Thus, with the above assumptions, the bulge may clearly 
be a cylindrical rotator in (N). In (K) the mean rotational 
velocity in the (£ = 5°, b = —20°) field is lower than that in 
the fields at (£ = 5°, 6 = —12°) and (£ = 5°, 6 = —15°) by 
Si 1.5 standard deviations, so cylindrical rotation cannot be 
ruled out. However, cylindrical rotation in (L) can be ruled 
out at the 2 a level. 

l.b.ii. With aRR : a 44 : Uzz in the ratio of solar neighbour¬ 
hood disk dispersions, 39 : 23 : 20, we find: 

(N) (£ = 5°, 6 = —12°): Qb = 20 ± 7kms“^kpc“^ and 
(7 = 72 ± 5 km s~^, 

(f=5°,6=—15°): Qb = 16±6kms~^kpc~^ and 
cr = 84 ± 10km s~\ 

(£ = 5° ,b = —20°): Qb = —7 ± 12 km s“^kpc“^ and 
(7 = 59 ± 7km s~^. 

(K) (f = 5°,6 = —12°): Qb = 30±llkms“^kpc“^ and 

(7 = 64 ± 5 km s“^, 

(f=5°,6=—15°): Qb = 15±8kms~^kpc~^ and 
(7 = 84 ± 7kms“\ 

(f=5°,6=—20°): Qb = 16±10kms~^kpc~^ and 
(7 = 82 ± 8 km s“^. 

(L) (f = 5°,6 = —12°): Qb = 20±7kms“^kpc“^ and 

(7 = 58 ± 6 km s~^, 

(f=5°,6=—15°): Qb = 14±9kms~^kpc~^ and 
(7 = 83 ± 8 km s“^, 

(f=5°,6=—20°): Qb = 16±9kms~^kpc~^ and 
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(7 = 81 ± 15 km . 

With the above assumptions, cylindrical rotation is consis¬ 
tent (at better than the 2a level) with scenarios (N), (K) 
and (L). 

l.b.iii. While if aRR : a^^ : Uzz is in the ratio of solar 
neighbourhood halo dispersions, 131 : 102 : 89, then 
(N)(l = 5°,6 = —12°): Qb = ll±5kms~^kpc~^ and 
(7 = 59 ± 3 km s“\ 

(1 = 5°, 6 = —15°): Qb = 16 ± 7km s~^kpc~^ and 
(7 = 77 ± 4 km s~^, 

(£= 5°, b = —20°): Qb = —7 ± 15 km s~^kpc~^ and 
(7 = 51 ± 6 km s“^. 

(K) (l = 5°,6 = —12°): Qb = 14±6kms~^kpc~^ and 

(7 = 61 ± 4 km s“\ 

(£ = 5°, b = —15°): Qb = 16 ± 9km s~^kpc~^ and 
(7 = 77 ± 5 km s~^, 

(f=5°,6=—20°): Qb = —7±14kms“^kpc“^ and 
(7 = 51 ± 7km s~^. 

(L) (1=5°,6 = —12°): Qb = 27±4kms“^kpc“^ and 

(7 = 58 ± 4 km s“\ 

(£ = 5°, b = —15°): Qb = 15 ± 10km s~^kpc~^ and 
(7 = 76 ± 6 km s“\ 

(£ = 5°, b = ~20°): Qb = —7±16 km s~^kpc~^ and 
(7 = 51 ± 6 km s“^. 

With the above assumptions, cylindrical rotation is consis¬ 
tent (at better than the 2a level) with scenarios (N), (K) 
and (L). 

According to the above hts to our data, cylindrical ro¬ 
tation is consistent with our data except if the bulge velocity 
ellipsoid is isotropic and the disk scale height is greater than 
165 pc. Cylindrical rotation in the bulge will be discussed 
in more detail below (where we will model the Galaxy less 
simplistically, dropping the assumption that the observed 
velocity distribution is only due to stars near the tangent 
point). 

In the paragraphs below, we will ht models to the ve¬ 
locity distributions of several helds at a time (see the intro¬ 
duction to this section). 

l.c. a isothermal, linear rotation curve 

We demand that the bulge velocity dispersion a is not a 
function of position in the Galaxy, he. da/dR = da/dz = 0 
and that it has linear rotation v{R) = £lbR (where Qb is 
constant). We select data from only those helds at 6 = —12° 
so as to suppress vertical dependence in the model hts. 

l.c.i. Constraining a such that the velocity ellipsoid is ev¬ 
erywhere isotropic, he. aRR = a,f„f, = azz, we hnd: 

(N) Qb = 27 ± 3 km s“^kpc“^ and cr = 58 ± 3 km s“\ 

(K) Qb = 27 ± 4 km s~^kpc~^ and cr = 57 ± 4 km s~^, 

This ht can be rejected (just) at the 2a level. 

(L) Qb = 25 ± 3 km s~^kpc~^ and cr = 56 ± 3 km s~^. 

This ht can be rejected at the 2a level. 

l.c.ii. Constraining aRR : a^^ : azz in the ratio of solar 
neighbourhood disk dispersions, 39 : 23 : 20, gives: 

(N) Qb = 27 ± 3 km s~^kpc~^ and aRR = 75 ± 5 km s~^, 

(K) Qb = 27 ± 2 km s“^kpc“^ and aRR = 74 ± 3 km s“\ 

(L) Qb = 25 ± 3km s“^kpc“^ and aRR = 73 ± 4 km s“^. 


All three hts can be rejected at the 2a level. 

l.c.iii. Next we try constraining aRR : a^^ : azz in the ratio 
of solar neighbourhood halo dispersions 131 : 102 : 89. This 
gives: 

(N) Qb = 27 ± 3 km s“^kpc“^ and aRR = 67 ± 4 km s“^, 
This ht can be rejected at the 2a level. 

(K) Qb = 27 ± 3 km s~^kpc~^ and aRR = 66 ± 4 km s~^, 
This ht can be rejected (just) at the 2a level. 

(L) Qb = 25 ± 2 km s~^kpc~^ and aRR = 65 ± 3 km s~^. 
Thus our data are consistent with the bulge being a lin¬ 
ear rotator in all the models considered above. The bulge 
rotation rate is well determined in all models. The velocity 
dispersion in these helds is much lower than that observed 
in the central regions of the bulge near Baade’s window, 
where the velocity dispersion of most tracer objects is in 
excess of lOOkms”^ (see e.g. the data compiled by Tyson 
1992). Tyson’s compilation of minor axis observations shows 
a clear decrease in velocity dispersion with increasing Galac- 
tocentric distance, with aRR tailing oh to 60 - 70kms~^ at 
|b| ^ 10. This is consistent with the above hts. The veloc¬ 
ity dispersion model which has the velocity ellipsoid in the 
shape of the solar neighbourhood disk dispersions is not con¬ 
sistent with the data (at the Si 2a level) — though it should 
be noted that the other hts are only marginally better (see 
appendix A). 

l.d. a isothermal, cylindrical rotation 

Next, we demand that the bulge velocity dispersion a is 
not a function of position in the Galaxy, i.e. da/dR = 
da/dz = 0 and also that the bulge behaves as a cylindrical 
rotator, whose mean rotation increases linearly with R over 
the region for which we have data, i.e. ir/{R,z) = UtR- 

l.d.i. Gonstraining a such that the velocity ellipsoid is ev¬ 
erywhere isotropic, i.e. aRR = a,f„f, = azz, we hnd: 

(N) Qb = 27 ± 2 km s“^kpc“^ and cr = 58 ± 3 km s“\ 

(K) Qb = 26 ± 2 km s“^kpc“^ and cr = 58 ± 4 km s“\ 

This ht can be rejected (just) at the 2a level. 

(L) Qb = 24 ± 2 km s“^kpc“^ and cr = 57 ± 3 km s“^. 

l.d.ii. Gonstraining aRR : a,f„f, : azz in the ratio of solar 
neighbourhood disk dispersions, 39 : 23 : 20, gives: 

(N) Qb = 27 ± 3 km s“^kpc“^ and aRR = 75 ± 2 km s“^, 

(K) Qb = 26 ± 2 km s“^kpc“^ and aRR = 76 ± 2 km s“^, 

(L) Qb = 24 ± 2 km s“^kpc“^ and aRR = 76 ± 4 km s“^. 

All three hts can be rejected at the 2a level. 

l.d.iii. Next we try constraining aRR : a,f„f, : azz in the 
ratio of solar neighbourhood halo dispersions 131 : 102 : 89. 
This gives: 

(N) Qb = 26 ± 2 km s“^kpc“^ and aRR = 67 ± 3 km s“^, 
This ht can be rejected (just) at the 2a level. 

(K) Qb = 26 ± 2 km s~^kpc~^ and aRR = 67 ± 2 km s~^, 
This ht can be rejected (just) at the 2a level. 

(L) Qb = 24 ± 2 km s“^kpc“^ and aRR = 67 ± 2 km s“^. 

It is clear from the above that the assumption of cylin¬ 
drical rotation, does not signihcantly ahect the htted rota¬ 
tion and velocity dispersion values listed in sub-paragraph 
(l.c) above. Again the velocity dispersion model which has 
the bulge velocity ellipsoid in the shape of the solar neigh¬ 
bourhood disk dispersions is not consistent with the data (at 
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the 2(7 level). The probabilities of the hts are summarised in 
Table 1. Thus, our data are (just) consistent with a model 
in which the bulge velocity dispersion is either isotropic or 
as anisotropic as the solar neighbourhood halo, and where 
the bulge is a rigid, cylindrical rotator. 

l.e. (7RR exponential, linear rotation 

Next we ask the question: assuming that the bulge is a 
linear rotator = ilbR, can the velocity dispersion 

be ht by an exponential function ^RR — ^Rr[R — 0,Z — 
f.Gkpc] exp( —(The scale length Rs is a free param¬ 
eter in the ht). This behavior of the velocity dispersion with 
R is observed in thin disk stars (see e.g. Gilmore, Wyse & 
Kuijken 1989). So as to eliminate most of the dependence 
of (7 on 2 in the hts, we consider only the four helds at 
h = -12°. 

l.e.i. We constrain a such that at a point, ctrr = = 

CTzz- The resulting hts have: 

(N) Qb = 24 ± 2 km s“^kpc“^, Rs = 5.2 ± 2.1 kpc, and 
itrrIR = 0,z = 1.6 kpc] = 71 ± 10 km s“\ 

(K) Qb = 24 ± 3 km s“^kpc“^, Rs = 4.6 ±1.8 kpc, and 
itrrIR = 0,z = 1.6 kpc] = 72 ± 8 km s“\ 

(L) Qb = 24 ± 3km s“^kpc“^, Rs = 3.3 ± 1.3kpc, and 
(TrrIR = 0, z = 1.6 kpc] = 77 ± 11 km s~^. 

l.e.ii. We also try a such that at a point, ctrr : : Uzz 

is in the ratio 39 : 23 : 20. Then: 

(N) Qb = 24 ± 2kms“^kpc“^, Rs = 22.5 ± 4.9kpc, and 
itrrIR = 0,z = 1.6 kpc] = 78 ± 9 km s“\ 

(K) Qb = 24 ± 3 km s“^kpc“^, = 18.2 ± 4.2 kpc, and 

itrrIR = 0,z = 1.6 kpc] = 79 ± 8 km s“\ 

(L) Qb = 24 ± 2 km s“^kpc“^, Rs = 8.6 ±2.7kpc, and 
itrrIR = 0,z = 1.6 kpc] = 82 ± 7 km s“^. 

This ht can be rejected at the 2(7 level. 

l.e.iii. We constrain a such that at a point, ctrr : : 

(Tzz is in the ratio 131 : 102 : 89. Then: 

(N) Qb = 24 ± 2 km s“^kpc“^, Rs = 8.2 ± 3.2 kpc, and 
itrrIR = 0, z = 1.6 kpc] = 76 ± 5 km s~^, 

(K) Qb = 24 ± 3 km s“^kpc“^, Rs = 7.0 ±2.3 kpc, and 
itrrIR = 0,z = 1.6 kpc] = 77 ± 10 km s“\ 

(L) Qb = 24 ± 3km s“^kpc“^, = 4.3 ± 1.5 kpc, and 

itrrIR = 0,z = 1.6 kpc] = 83 ± 8 km s“^. 

It is interesting to note that, especially in the isotropic 
velocity dispersion situation, the htted radial scale length Rs 
is broadly similar to that of the thin disk and that assumed 
for the thick disk. [Recall that simple dynamical models for 
the disk suggest that the scale length for variation of the 
stellar velocity dispersion is just twice that for variation of 
the spatial density.] Thus with the assumption that <trr(R) 
is an exponential function, it is possible to ht the bulge with 
velocity dispersion models where the velocity ellipsoid is as 
anisotropic as that of the solar neighbourhood disk. 

In summary, our data set implies a well constrained lin¬ 
ear rotation curve in the bulge, but it allows many bulge ve¬ 
locity dispersion relations; the shape of the velocity ellipsoid 
can be isotropic or as battened as the solar neighbourhood 
disk, and there may, or may not be, radial variations in the 
velocity dispersion. This is discussed further in section 11. 


2. Can the outer bulge be triaxial? 

Though many studies have found that the bulge is non- 
axisymmetric (c/ section 1), it should be stressed that the 
data we have obtained can be adequately described by sim¬ 
ple axisymmetric models. Perhaps when the velocities of 
many more ‘outer bulge’ K giants are known, it may be 
necessary to resort to a less simple explanation. Neverthe¬ 
less, we investigate whether our data are compatible with 
non-axisymmetric models. 

Published non-axisymmetric bulge models vary signifi¬ 
cantly according to the bulge tracer they are fit to, and so 
are not well established. To reduce the large numbers of free 
parameters, we assume the the bulge has a triaxial form sim¬ 
ilar to that found by Binney et al.(1991) (which was fit to the 
central Galactic bar): the axial ratios are 6/a = c/a = 0.75 
and the bar is viewed at an angle of 16° to its major axis. 
We adopt the density profile used by Binney et al.(1991): 
ic(s) = r(o(s/so)~^'^, where ± (y/0.75)^ ± (2/0.75)^ 

and So = 1200 pc. We further assume that the closed ‘ri’ 
orbits (Gontopoulos & Mertzanides 1977) follow approxi¬ 
mately the isodensity surface, and that the velocity ellipsoid 
is aligned with it such that one diagonalising axis (u) points 
out normal to the isodensity surface, another diagonalising 
axis (c) points along the isodensity surface in a plane paral¬ 
lel to the Galactic plane, while the third diagonalising axis 
(re) is perpendicular to the other two. We try models where 
bulge rotation is a linear function of Galactocentric radius 
R (independent of ^ and bulge geometry). 

2. a. With the above assumptions we can reject at the 4(7 
level the most likely fit to a model that has v,f, = QbR, 
^uu — ^vv — ^ww and (7^i^i(s) — constant. 

2.b. Assuming also that v,f, = QbR, (7„„(s) = constant, 
(7u,uj(s) = constant and = constant, we hnd: 

(N)Qb = 25.4± 4.1 km s“^kpc“^, 

(7im = 87 ± 9 km s“^, cTzz = 30 ± 6 km s“^ and (7u,u, = 
50 ± 11 km s“\ 

(K) Qb = 25.0±3.9 km s“^kpc“^, 

(7im = 87 ± 7km s“^, cTzz = 28 ± 6 km s“^ and (7u,u, = 
50 ± 12km s~\ 

(L) Qb = 24.8±3.5kms“^kpc“^, 

(7im = 81 km s“^ 14, (7 1 , 1 , = 24±7kms“^ and (7u,u, = 
84 ±48 km s~P 

These hts require a very non-isotropic velocity ellipsoid. 
This would imply a signiheant change in the shape of the 
velocity ellipsoid from our innermost held to Baade’s win¬ 
dow, where the velocity dispersion is observed to be isotropic 
(Spaenhauer et al.l992). This calls the physical plausibility 
of this ht into question, unless the ‘nuclear bulge’ is unre¬ 
lated to the ‘outer bulge’. 

2.C. We can reject, again at the 4(7 level, the most likely 
ht to a model that has v,f, = QbR, <ruu = <rvv = <rww and 
^uuG) = ^uu[^ = 0] exp( —s/2ss) (where Ss is a variable). 

2.d. An isotropic velocity ellipsoid can be ht to our data, if 
we free the orientation a of the prolate bulge: 

(N)Qb = 24±4.2 km s“^kpc“^, 

(7 = 56 ± 7 km s~^ and a = 29° ± 6°, 

(K)Qb = 24±2.7 km s“^kpc“^, 

(7 = 54 ± 5 km s~^ and a = 26° ± 5°, 
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(L) Qb = 25.2 ± 3.2 km s“^kpc“^, 

(7 = 60 ± 5 km and a = 27° ± 4°. 

Though a = 27° ± 5° is at odds with the central bar 
model of Binney et al.(1991), it is consistent with the triaxial 
HI model of Blitz & Spergel (1991). Thus, with several 
simplifying assumptions, it is possible to ht our data to bulge 
models with non-axisymmetric geometry. 

The last example above demonstrates the problem of 
htting non-axisymmetric models to our data: the constraints 
on the geometry of the bulge from extant studies allow too 
great a range of models. 

8 SPECIFIC ANGULAR MOMENTUM 
DISTRIBUTION OF THE BULGE 

The distribution of specihc angular momentum is a primary 
attribute of a Galactic component (see Wyse & Gilmore 
1992). To obtain this distribution one clearly needs the 
rotation and density distribution through the component. 
We established in the previous section that the rotation 
in the bulge from 700 pc < R < 3500 pc at 6 = —12° is 
very well constrained. Almost independent of the model 
assumed for the velocity dispersion, the bulge rotates at 
Qb = 25 km s~^kpc~^. Kent (1992) gives a spatial density 
model of the bulge which is ht to 2.4jUm emission for the 
outer bulge regions and to the data of Allen et al.(1983) 
and Becklin & Neugebauer (1968) for the central regions. 
The model has luminosity density n, such that: 

n = 1.04 X 10® (s/0.482)“^ “ Lopc"® 

for s < 938 pc and 

n = AT (s/667) Lqpc”® 

for s < 938 pc. In this ht s is an oblate cylindrical coordinate 
such (s^ = -f ( 2 / 0 . 6 !)^) and AT is a modihed Bessel 
function. 

This density distribution is integrated through the 
Galaxy, to give the mass M out to radius R. From M{R) one 
can directly obtain the distribution of specihc angular mo¬ 
mentum M{h) by substituting for the angular momentum 
h = HbR^ ■ The distribution of specihc angular momentum 
of the bulge is displayed in Figure 8 , where we also show the 
specihc angular momentum distribution of the halo, disk 
and thick disk, as given by Wyse & Gilmore (1992). The 
bulge, under the above assumptions, has an even more dis¬ 
sipated specihc angular momentum distribution than found 
by Wyse & Gilmore (1992). This is of course consistent with 
the hypothesis that the bulge is the dissipated core of the 
halo, but is not consistent with the hypothesis that the bulge 
is closely related, in an evolutionary sense, to the thick disk 
or the thin disk. 


9 THE BULGE METALLICITY 
DISTRIBUTION FUNCTION 

We show a representative colour-magnitude (GM) diagram 
of the helds studied. Figure 9 is the GM diagram in (Bj — R) 
and R as found for the (£ = 5°, 6 = —12°) held. The ‘star’ 
symbols represent the stars observed spectroscopically. In 
this held, the distance modulus to the tangent point to 


the line of sight is (m — M) = 14.5, and the reddening is 
Eb-v = 0.14. By superimposing GM relations for four 
globular clusters of diherent metallicities (shifted to the tan¬ 
gent point), we illustrate the metallicity range of stars that 
have been observed near the bulge. (The colour-magnitude 
system (Bj — R), R was converted into (B — V), V using 
the colour equations presented in Paper 1). The clusters are 
from left to right: M92, M3, M5, 47 Tucanae and NGG 188. 
At the tangent point in this held we therefore expect to 
hud few halo stars signihcantly more metal poor than M3 
([Fe/H] = —1.66 - Djorgovski 1993), but more metal rich 
populations should easily be detected, if present. 

In a high-resolution spectroscopic study of bulge K- 
giants, McWilliam & Rich (1994) hud that Mg and Ti in 
their stars are more abundant by Si 0.3 dex than in the Sun 
over the range of [Fe/H] they observed, while Ga and Si 
are present in similar ratios to those of disk giants. If that 
deduction is correct, we cannot simply infer [Fe/H] from 
[Mg/H], as measured by the Mg index dehned in Paper 1, 
which was calibrated against local K-giants. Further data 
would be necessary to establish the bulge [Mg/Fe] gradient. 
When the correctness of the bulge element ratio results has 
been established, a correction will need to be applied to the 
abundances derived here. 

With the above health warning, we show the relation¬ 
ship between metallicity and radial velocity for objects clas- 
sihed as K giants in Figure 10. Recall that reliable cal¬ 
ibration of the photometry is not available for the helds 
(i = —15,6 = —12) and (f = 5,6 = —20), so that for these 
helds [Fe/H] is really only a ranking in order of relative 
abundance. Furthermore, we show in IGl that the low ve¬ 
locity dispersion feature with mean velocity x 170kms~^ 
seen in Figure lOd, Figure lOe and Figure lOf is the Sagit¬ 
tarius dwarf, an object that is not part of the Milky Way, 
at least, just yet. We therefore can conhdently analyse the 
abundance - radial velocity data of Figure 10a, Figure 10c, 
Figure lOd and Figure lOe , and for the latter two helds we 
select only those stars with Galactocentric radial velocity 
V < 130kms~^: we proceed to do this in the remainder of 
this section. 

Recall that Si 90% of objects in our sample classihed as 
K dwarfs have formally-calculated [Fe/H] > 0.0. The tech¬ 
nique for determining these values was calibrated on K gi¬ 
ants only, so for K dwarfs the calculated ‘metallicity’ val¬ 
ues will be very diherent from the actual metallicity values. 
In most helds, the expected number of dwarf stars will be 
Si 20%, while (mainly) photometric errors will induce Si 25% 
misclassihcation. Given that Si 20% of objects classihed 
as K giants have [Fe/H] > 0.0, we should expect Si 25% K 
dwarf contamination in the K giant sample for stars with 
calculated [Fe/H] > 0.0. Thus, the observed metallicity dis¬ 
tributions are unreliable for [Fe/H] > 0.0. 

The velocity dispersion in the four helds at {£ = 
-25°,6 = -12°), (£ = -5 °,6 = -12°), (£ = 5 °,6 = -12°) 
and (f = 5°, 6 = —15°) decreases with increasing metallicity. 
We divide the stars into three bins; these are dehned as: 

i) bin-1. [Fe/H] < —1.68 — that is, a region containing 
stars more metal poor than the mean halo abundance 
(Laird et al.l988) 

ii) bin-2. 1.68 < [Fe/H] < 0.0 — the region between the 
mean halo abundance and the abundance where we ex¬ 
pect dwarf and distant disk giant contamination to be- 


16 R. Ibata and G. Gilmore 


Table 1: KS test probabilities for the bulge functional fits for all the models described in the text. 


Test 

(-25°, -12°) 

(-15°, -12°) 

(-5°, -12°) 

(5°,-12°) 

(5°, -15°) 

(5°, -20' 

l.a.i 

0.52 

0.69 

0.47 

— 

— 

— 

l.a.ii 

0.45 

0.53 

0.58 

— 

— 

— 

l.a.iii 

0.46 

0.41 

0.66 

— 

— 

— 

l.b.i(N) 

— 

— 

— 

0.43 

0.11 

0.96 

l.b.i(K) 

— 

— 

— 

0.45 

0.13 

0.97 

l.b.i(L) 

— 

— 

— 

0.33 

0.16 

0.96 

l.b.ii(N) 

— 

— 

— 

0.43 

0.14 

0.96 

l.b.ii(K) 

— 

— 

— 

0.12 

0.15 

0.14 

l.b.ii(L) 

— 

— 

— 

0.55 

0.21 

0.17 

l.b.iii(N) 

— 

— 

— 

0.38 

0.14 

0.97 

l.b.iii(K) 

— 

— 

— 

0.12 

0.15 

0.97 

l.b.iii(L) 

— 

— 

— 

0.23 

0.11 

0.97 

l.c.i(N) 

0.13 

0.28 

0.10 

0.29 

— 

— 

l.c.i(K) 

0.11 

0.25 

0.07 

0.35 

— 

— 

l.c.i(L) 

0.08 

0.11 

0.03 

0.40 

— 

— 

l.c.ii(N) 

0.28 

0.41 

0.07 

0.43 

— 

— 

l.c.ii(K) 

0.26 

0.37 

0.08 

0.22 

— 

— 

l.c.ii(L) 

0.17 

0.21 

0.04 

0.27 

— 

— 

l.c.iii(N) 

0.19 

0.34 

0.05 

0.22 

— 

— 

l.c.iii(K) 

0.16 

0.31 

0.07 

0.27 

— 

— 

l.c.iii(L) 

0.10 

0.15 

0.18 

0.33 

— 

— 

l.d.i(N) 

0.43 

0.19 

0.28 

0.42 

0.15 

0.12 

l.d.i(K) 

0.20 

0.27 

0.09 

0.32 

0.08 

0.09 

l.d.i(L) 

0.10 

0.12 

0.15 

0.37 

0.10 

0.11 

l.d.ii(N) 

0.46 

0.41 

0.06 

0.15 

0.08 

0.07 

l.d.ii(K) 

0.41 

0.38 

0.12 

0.20 

0.09 

0.05 

l.d.ii(L) 

0.27 

0.21 

0.07 

0.24 

0.11 

0.04 

l.d.iii(N) 

0.32 

0.35 

0.10 

0.20 

0.07 

0.15 

l.d.iii(K) 

0.28 

0.32 

0.13 

0.26 

0.10 

0.08 

l.d.iii(L) 

0.16 

0.16 

0.18 

0.30 

0.11 

0.19 

l.e.i(N) 

0.97 

0.13 

0.18 

0.49 

— 

— 

l.e.i(K) 

0.99 

0.22 

0.15 

0.53 

— 

— 

l.e.i(L) 

0.97 

0.34 

0.06 

0.67 

— 

— 

l.e.ii(N) 

0.71 

0.12 

0.20 

0.71 

— 

— 

l.e.ii(K) 

0.92 

0.25 

0.17 

0.45 

— 

— 

l.e.ii(L) 

0.96 

0.29 

0.11 

0.37 

— 

— 

l.e.iii(N) 

0.89 

0.11 

0.20 

0.43 

— 

— 

l.e.iii(K) 

0.99 

0.23 

0.16 

0.47 

— 

— 

l.e.iii(L) 

0.98 

0.29 

0.18 

0.42 

— 

— 

2.b(N) 

0.87 

0.36 

0.16 

0.66 

0.34 

0.33 

2.b(K) 

0.65 

0.41 

0.19 

0.78 

0.35 

0.29 

2.b(L) 

0.51 

0.20 

0.11 

0.85 

0.38 

0.31 

2.d(N) 

0.87 

0.36 

0.16 

0.66 

0.34 

0.33 

2.d(K) 

0.65 

0.41 

0.19 

0.78 

0.35 

0.29 

2.d(L) 

0.42 

0.15 

0.18 

0.71 

0.23 

0.49 
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Figure 8. The specific angular momentum distributions for the four major Galactic components. The solid curve is the distribution of 
the bulge component, as given by the rotational velocity found in this work together with a density model of Kent (1992). The other 
curves are taken from Wyse &: Gilmore (1992); the dashed-dotted curve represents the halo, the dotted curve represents the thick disk 
and the dashed curve represents the disk. 


Figure 9. (Bj — R, R) colour-magnitude diagram of stars mea¬ 
sured in a 1° X 1° region of sky in the [i = 5°, b = —12°) field. 
The ‘star’ symbols represent the stars observed spectroscopically. 
The four colour-magnitude relations superimposed are, from left 
to right: M92, M3, M5, 47 Tucanae (from Sandage 1982) and 
M67 (from Chiu 1980), all shifted to the distance of the tangent 
point in this line of sight. 

come significant 

hi) bin-3. [Fe/H] > 0.0 — where we expect significant 
dwarf and distant disk giant contamination 


The velocity dispersions in Figure 10a are found to be 
85.7kms“\ 59.5kms“^ and 46.4kms“^ in the first, second 
and third bins respectively, while in Figure 10c the velocity 
dispersions are 112.3 km s~^, 70.9 km s~^ and 58.3 km s~^ in 
the hrst, second and third bins respectively. We do not 
calculate similar quantities in Figure lOd and Figure lOe, 
because of the contamination in those helds. 

We model the metallicity distribution of the halo 
as a gaussian of intrinsic width a = 0.2dex centred at 
Fe/H] = —1.68 (Laird et al.l988), and we assume that the 
shape of the thick disk metallicity distribution towards the 
South Galactic Pole (Gilmore, Wyse & Jones 1995) is the 
same as that in the helds studied here. That is, we assume 
that there are no radial abundance gradients in either the 
halo or the thick disk. So we pose the question: is it possible 
to ht the observed metallicity distributions with a composite 
model that consists of the sum of the above halo and thick 
disk models together with a ‘simple model’ of chemical evo¬ 
lution to represent the bulge? 

We hrst need to know the relative numbers of halo and 
thick disk stars in our sample: in table 2, we show the num¬ 
ber of halo and thick disk stars predicted by the galaxy 
model in the four helds investigated. 

The simple, or closed-box, model of Galactic chemical 
evolution is discussed in standard texts, and is discussed by 
Rich (1990) in the context of the Milky Way bulge. The 
probability density function dW/dpe/H] it predicts is given 
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Figure 10. The relationship between velocity and metallicity for stars classified as giants. The fields are: (a) i = —25°, h = —12°, (b) 
I = -15°, b = -12°, (c) i = -5°, b = -12°, (d) ^ = 5°, 6 = -12°, (e) ^ = 5°, 6 = -15° and (f) ^ = 5°, 6 = -20°. 


by 


dN 

rf[Fe/H] 


y y 


(6) 


where y is the yield and zq is the solar metal abundance. In 
an evolved population which has turned its gas into stars, 
the mean abundance equals the yield, i.e. 1 = y. 

Since we are not conhdent of the observed distributions 
in the region [Fe/H] > 0.0, we will truncate both the model 
and the data at [Fe/H] = 0.0. We stress that this cutoff 
will have the effect of altering the yield of the resulting hts. 


However, it is useful to have a simple parameterisation of the 
data, and this specihc limitation in application is relatively 
minor compared to the rather extreme nature of some of the 
physical assumptions underlying the simple model. 

To simulate measuring errors in the above composite 
model, we convolve the empirical thick disk distribution of 
Gilmore, Wyse & Jones (1995) with a gaussian of cr = 0.25 
(the difference of the measuring errors between their and 
our [Fe/H] measurements), while the theoretical bulge and 
halo distributions are convolved with our metallicity error 
distribution (c/ Paper 1). The only free parameter is the 
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Table 2. The predicted number of halo and thick disk stars with 
[Fe/H] < 0.0 according to the Galaxy model. In the (5°,—12°) 
and (5°, —15°) fields, we have also only selected stars with Galac- 
tocentric radial velocity v < 130kms“^. 


Component 

(-25°, -12°) 

(-5°,-12°) 

(5°,-12°) 

(5°,-15°) 

Total 

158 

206 

219 

66 

Halo 

3.5 

4.3 

4.8 

1.2 

Thick disk 

32.2 

30.3 

32.2 

5.6 


yield of the ‘simple model’, which we hnd by maximising 
the likelihood of the ht. The results are shown in Figure 11; 
all these hts are acceptable at better than the l.Scr level. The 
yield y in units of zq is found to be 0.69 ± 0.20, 0.43 ± 0.11, 
0.53 ±0.06 and 0.46 ±0.09 in Figure 11a, Figure 11b, Figure 
11c and Figure lid respectively. The reader can use the 
agreement or otherwise between these derived yields as an 
indication of the ability of the model to represent the data 
adequately. 

The kinematic information can now be used to check for 
self-consistency in this model. As usual, we will not use the 
helds at (f = 5°, 6 = —12°) and (£ = 5°, b = —15°), because 
of the contamination from the dwarf galaxy. Note that ac¬ 
cording to the adopted thick disk metallicity distribution, 
less than 0.5% of stars in bin-1 are thick disk stars, while 
thin disk giants comprise less than 5% of the stars in bin-2 
and contribute negligibly to bin-1. We then ht the kinemat¬ 
ics of bulge stars in bins 1 and 2 separately using our Galaxy 
model; standard kinematic and space density distributions 
are assumed for the halo and thick disk (with normalisations 
in bins 1 and 2 as detailed above), while the bulge compo¬ 
nent is modeled as an oblate rotator with isotopic velocity 
dispersion throughout. Descriptions of these distributions 
may be found in Gilmore, King & van der Kruit (1990) and 
in Gilmore, Wyse & Kuijken (1989). The bulge is assumed 
to rotate rigidly at a rate Q = 25kms~^ kpc~^ (as derived 
in §7.3), which leaves the bulge dispersion as the only free 
parameter in the ht. 

A maximum likelihood ht to the data in the {£ = 
— 25, h = —12) held gives a = 82±46 in bin-1 and a = 69±7 
in bin-2, while a ht to the data in the {£ = —5,6 = —12) 
held gives a = 139 ±34 in bin-1 and cr = 53±6 in bin-2. We 
therefore conclude that either there are signihcantly more 
halo stars in the held at {£ = —5,6 = —12) than the stan¬ 
dard halo model predicts (so as to increase the ratio of halo 
to bulge stars in bin-1), or there is a signihcant dependence 
of velocity dispersion on metallicity in the bulge. 

Alternatively, if we demand that there is no change in 
bulge dispersion with metallicity in the {£ = —5,6 = —12) 
held, then there must be approximately 25 halo stars in the 
sample (a factor of Si 5 more than the standard model), and 
the bulge dispersion becomes 99±30kms“^, consistent with 
53 ± 6kms~^. If this assumption is true, we hnd that the 
composite metallicity model described above, which has a 
‘simple model’ with yield y = 0.72 ± 0.26 representing the 
bulge component, is acceptable at better than the la level; 
the ht is shown in Figure 12. 

How much are the above hts ahected by our decision 
to ignore bin-3 data? We take the numbers of halo and 


thick disk stars as given in table 2, then ht the composite 
model to the data in all 3 bins. The resulting yield (in units 
of Zq) is found to be y = 0.76 ± 0.07, y = 0.40 ± 0.05, 
y = 0.64 ± 0.08 and y = 0.74 ± 0.10 in the held at (£ = 
-25°,b = -12°), (£ = -5°,6 = -12°), (£ = 5°,6 = -12°) 
and {£ = 5°,b = —15°) respectively (the distributions are 
acceptable at better than the l.5a level). Gomparison with 
the results above shows that imposing the abundance cut 
marginally reduces the derived yield, as expected. 

In summary, the yield found in each held is dependent 
on the imposed metallicity cut and slightly dependent on the 
assumed normalisation of thick disk and halo stars. Never¬ 
theless, it appears that the bulge component can be ade¬ 
quately well described (for [Fe/H] < 0.0) by the abundance 
distribution predicted by the ‘simple’ or closed box model 
of chemical evolution. From htting this simple model we 
obtain the yield of a simple-model ht to the K giant popula¬ 
tion in the neighbourhood of the tangent point in each held: 
this is consistent with no metallicity gradient through the 
region in which we have data: i.e. 0.7kpc D, R Iz 3.5 kpc. 
Recall also that recent developments (McWilliam & Rich 
1994) have shown that K giants in Baade’s window typi¬ 
cally have abundance [Fe/H] Si —0.25 — this would imply 
a zero metallicity gradient throughout almost the whole of 
the bulge. 


9.1 A Connection with the Thick Disk? 

Since the thick disk metallicity distribution in Figure 11 
looks similar to that of the bulge, we will investigate whether 
we can ht the observed metallicity distributions with a 
model similar to that discussed above, but whose bulge com¬ 
ponent has the same metallicity distribution as the thick 
disk. 

We hnd a good (better than 1.5(7) ht by adopting the 
Galaxy model halo star prediction of 4 halo stars in the 
(£ = —25°, 6 = —12°) held; this is displayed in Figure 13a. 
If the model predictions for the number of halo stars in the 
other helds are correct, then we can reject at the 2a level the 
hypothesis that the bulge has the same abundance distribu¬ 
tion as the thick disk in those helds. To obtain agreement 
at better than the 1.5(7 level, we require 30 halo stars in the 
(£ = —5, 6 = —12) held, 40 halo stars in the (7=5,6= —12) 
held and 20 halo stars in the (7 = 5,6 = —15) held. Recall 
that 4 such stars are expected. These hts are shown in Fig¬ 
ure 13. 
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Figure 11. The metallicity distributions in (a) the field at (t = —25°, 6 = —12°), in (b) at (t = —5°, 6 = —12°), in (c) at 
[i = 5°, 6 = —12°) and in (d) at [i = 5°, 6 = —15°). These are compared to the model described in the text, in which the bulge 
component is represented by a simple ‘closed box’ chemical evolution model, with a yield fitted independently to each field. In each 
panel, the dashed line represents the halo component, the dashed-dotted line represents the thick disk, the dotted line represents the 
bulge and the full line is the sum of all three components. 


10 COMPARISON TO EXTERNAL GALAXIES 

Kormendy & Illingworth (1982; hereafter KI82) compared 
the importance of rotational support in spiral galaxy bulges 
to that in elliptical galaxies. The ratio of the energy in or¬ 
dered motion to that in random motion was estimated from 
the ratio of the maximum rotational velocity to the cen¬ 
tral velocity dispersion: Vm/cro. Vm/cro is approximately 
equal to the mass-weighted averaged quantities {V^ 
in isotropic oblate spheroid models (Binney 1980). It is the 
mean stellar rotation that provides support and flattening 
in these models. Similarly, Davies et al.(1983a) showed that 
the kinematics of elliptical galaxies of low absolute luminos¬ 
ity are also well described by such models, while intrinsically 
bright elliptical galaxies bear most resemblance to models 
which are supported by random motion (ie pressure) and 


flattened by anisotropic velocity dispersions. 

An interesting feature of the data of KI82 is the lack of 
spatial overlap between their data for external bulges with 
extant data for the Milky Way. The only regions of the 
Milky Way not dominated by the disk where the kinemat¬ 
ics of K-giants were well known are the Baade’s Window 
bulge and the subdwarf halo in the solar neighbourhood. 
Comparable regions in the KI82 edge-on spiral galaxies were 
either hidden by the galactic disks or were too faint to be 
detected. KI82 quote the limiting surface brightness in their 
data as 24 B mag arcsec”^, while Morrison (1993) estimates 
that the surface brightness of the halo, as viewed edge-on 
from outside the Galaxy at the solar radius, is jUv = 27.7 V 
mag arcsec”^. The various contributions to the light profile 
observed by KI82, and the consequences thereof, were dis¬ 
cussed in some detail by Shaw & Gilmore (1989). The essen- 
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Figure 12. Comparison of the data towards (t = —5°,6 = —12°) 
with the composite model described in the text. We have as¬ 
sumed that there are 25 halo stars in the sample rather than 
the 4 expected. The dashed line represents the halo component, 
the dashed-dotted line represents the thick disk, the dotted line 
represents the bulge and the full line is the sum of all three com¬ 
ponents. 

tial point is that no data were available to compare available 
spectroscopic studies of external spiral bulges with the bulge 
of the Milky way. We now have the necessary kinematic data 
to compare the Milky Way bulge to the spiral galaxy bulges 
observed by K182, over a similar range of R/Hr and zjhz, 
where tiR and hz are characteristic bulge scale lengths in the 
R, z directions. We also investigate whether the Milky Way 
bulge is typical of spiral galaxy bulges. Since the bulge spa¬ 
tial density prohle is now reasonably well known (e.g. Kent 
1992; Freeman 1993), we also investigate whether the above- 
mentioned assumption, used by K182, that Vm/cro ~ Km/^, 
holds true in our Galaxy. 

In most of the galaxies they studied, K182 were not 
able to measure Vm and ao directly due to problems with 
disk contamination, but used instead interpolations to es¬ 
timate those parameters. We reproduce their technique on 
the Milky Way bulge. 

K182 calculated Vm from those bulge regions where 
“V (r) reaches a well-dehned plateau or maximum at radii 
dominated by bulge light”. This is the case in our (I = 
— 25°, 6 = —12°) held, where the mean line of sight ve¬ 
locity has increased by less than 15kms“^ over a range of 
Si 1 kpc from (I = —15°, 6 = —12°). Since the mean rota¬ 
tion in most galaxy bulges decreases with height above the 
galactic plane, K182 extrapolated the Vm values down to 
Vm{z = 0). The extrapolation was performed by using the 
functional form of V (z) as obtained from two-dimensional 
velocity maps of elliptical galaxies by Davies and Illingworth 
(1983b). V (z) from these elliptical galaxies was scaled by the 
length Ze, the value of the effective radius of a Vaucouleurs 

law ht to the minor axis prohle. We use Ze = 200pc for 
the Milky Way bulge from Frogel et al.(1990). Where the 
external galaxy bulge displayed cylindrical rotation, K182 
used Vm{z = 0) = Vm- However, whether or not the 
Milky Way bulge is a cylindrical rotator is not as yet estab¬ 
lished (c/section 4), so we will treat both cases separately. 
Cylindrical rotation would then give Vm{0) = 95kms“^ at 
{£ = —25°, h = 0°), while applying the above V(z) correction 


for the non-cylindrical rotation case, we obtain Vm{0) = 104. 

K182 approximate central dispersion ao to an average W 
out to Si re/2, where re is the ehective radius of a de Vau¬ 
couleurs (1959) law ht. No K giant data are available 
over the same Galactocentric distance range. However, cen¬ 
tral Miras r Si 350 pc from the Galactic centre have velocity 
dispersion 127 ±24kms“^ Feast et al.(1990). This value is 
compatible with the velocity dispersion of 131±28kms~^ of 
Carbon stars at r Si 450pc found by Tyson & Rich (1991). 

Thus Vm/cTo = 0.75 if the bulge is a cylindrical rotator, 
or Vm/cTo = 0.82 if it resembles the other model above. 

Therefore, if the bulge is similar to the oblate rotating 
models with isotropic velocity dispersion (Binney 1978), to 
which K182 compare their data (see K182, their hgure 3), 
then we hud that the ellipticity of the Milky Way bulge 
should be e Si 0.4. This value of bulge ellipticity is very 
similar to the value e = 0.39 found by Kent et al.(1991). 

11 CONCLUSIONS 

Our analysis of the data discussed in Paper 1 has shown that 
the outer regions of the Galactic bulge between 0.7kpc < 
R V 3.5kpc and 1.2 < \z\ < 2.7kpc can be adequately de¬ 
scribed by axisymmetric models with linear rotation Q = 
25 km s~^ kpc~^. The data also constrain the Galactic disk 
to have vertical scale height Zh ^ 220kms~^ in the central 
Si 3.5 kpc. Since the vertical scale height of the Galactic 
disk in the solar neighbourhood is Zh = 250 pc (Kuijken 
& Gilmore 1989), the Galactic disk appears to be different 
from pure-disk galaxies, where the disk scale height is inde¬ 
pendent of galactocentric distance. Alternatively, our selec¬ 
tion function for K giants continues to select preferentially 
in favour of the young disk in preference to the old disk. 
The data allow several velocity dispersion cr{R) models. If 
cr{R) is isothermal, then the bulge velocity ellipsoid can only 
be about as anisotropic as that of the solar neighbourhood 
halo. However, if cr{R) decreases exponentially with R, the 
bulge velocity ellipsoid can be as anisotropic as that of the 
solar neighbourhood disk. More data in several more helds 
will be required to conhrm or reject cylindrical rotation and 
to constrain better the radial dependence of the shape, size 
and orientation of the velocity ellipsoid. The simplest al¬ 
lowed model which is consistent with this data-set does not 
require non-axisymmetric geometry. 

We hnd acceptable agreement with the axisymmetric 
bulge model devised by Kent (1992), which uses a Galactic 
mass model ht to 2.4jUm luminosity, with the assumption 
M/L = 1, to predict self-consistently the kinematic struc¬ 
ture of the bulge. Thus the potential of the inner Galaxy is 
currently reasonably well established. 

We show that the shape of the observed abundance dis¬ 
tribution (for [Fe/H] < 0.0, where we believe are data are 
reliable) can be adequately represented by a model that con¬ 
sists of the sum of the abundance distributions for the halo 
and thick disk observed in the solar neighbourhood, together 
with an abundance distribution for the bulge as predicted 
by the ‘simple model’ of chemical evolution. The maximum- 
likelihood htted yield in each held does not vary signihcantly 
over the Galactic regions for which we have data, and is con¬ 
sistent with a zero metallicity gradient in the bulge K giant 
population from the Baade’s Window region to R K, 3.5 kpc. 
Thus the bulge is a well-mixed population. 
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Figure 13. We test the hypothesis that the bulge has a similar metallicity distribution to the empirical thick disk distribution at the 
Solar Galactocentric distance determined by Gilmore, Wyse & Jones (1995). The metallicity distributions correspond to (a) the field at 
(I = —25°, h = —12°), in (b) to (I = —5°, h = —12°) in (c) (£= 5°, b = —12°) and in (d) (£ = 5°, b = —15°). The number of halo stars 
required to make these fits acceptable at better than the 1.5(T level, are 3.5, 30, 40, and 5 stars in (a), (b), (c) and (d) respectively, while 
~ 4 are expected in each field. The dashed line represents the halo component, the dotted line represents the bulge and the full line is 
the sum of the two components. 


The bulge region investigated in this volume covers ap¬ 
proximately the same Galactocentric distance range (mea¬ 
sured in bulge scale lengths) as did the kinematic study 
of external spiral galaxy bulges presented by Kormendy & 
Illingworth (1982). Our data allow the first detailed kine¬ 
matic comparison of the Milky Way bulge to the population 
of external spiral galaxy bulges. We find that the Galactic 
bulge is representative of the population of external spiral 
galaxy bulges. 

With our measured linear rotation curve, and adopting 
the mass model derived from the observed near-lR surface 
brightness of the bulge by Kent (1992), we derive the specific 
angular momentum distribution of the bulge. This distribu¬ 
tion is very similar to, though slightly more dissipated than, 
that of the halo, and so is consistent with the hypothesis 


that the bulge is the dissipated core of the halo. The cor¬ 
responding distributions of specific angular momentum for 
the disks are very different. Thus, our results are not eas¬ 
ily consistent with any model in which the bulge is closely 
related, in an evolutionary sense, to the Galactic disks. 

Finally, we note that our study is of what we have 
termed the Outer Bulge. Although this includes most of the 
volume of the central Galactic bulge, in the sense that this 
term is applied to external galaxies, most recent studies of 
the central Galaxy have been more concerned with the very 
inner parts of the Bulge. This part of the Galaxy is nicely 
seen in recent IRAS and GOBE maps, but is nearly com¬ 
pletely obscured from short wavelength, star by star study. 
This very central, and very compact, structure is the central 
highest density part of all Galactic structural components. 
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It is possible that something else also lives there. We are 
conhdent that the order of magnitude improvements in spa¬ 
tial resolution and sensitivity which will be provided by ISO 
at wavelengths appropriate to study the dominant central 
stellar population will provide substantial improvements in 
our knowledge of this very central Galactic core. In the in¬ 
terim, we refrain from speculation as to the relationship, if 
any, between that region and the outer bulge which has been 
the subject of this investigation. 
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